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Abstract 

Based on a Lagrangian with a coupling of two gluons to J 71 " = + (the quan- 
tum numbers of the vacuum) which decay to qq pairs, a model is presented, in 

+- > ■ 

which hadrons couple directly to the absolute vacuum of fluctuating gluon fields. 

O : 

By self-consistency requirements the confinement potential as well as qq densities 
and masses are obtained, which are in good agreement with experimental data on 
4^3 | scalar and vector mesons. In comparison with potential models additional states 

or 

43 



(N 



are predicted, which can explain the large continuum of scalar mesons in the low 
mass spectrum and new states detected recently in the charm region. 
The presented model is consistent with the concept, that the hadron masses can be 



understood by binding effects of the quarks. 



PACS/ keywords: 12.38 Aw, 12.38.Lg, 12.39.Mk/ Gluon-gluon coupling to qq- 
pairs, generation of stable gg-systems assuming massless quarks, confinement po- 



\Q ■ tential, masses of ++ and 1 mesons. 

o 
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In the hierarchy of quantum systems hadrons represent the smallest complex substructures 
known inside of atoms and nuclei. This is supported by the property of asymptotic 
freedom [1] of the strong interaction. Therefore, hadrons can be related directly to the 
absolute vacuum of fluctuating gluon fields (with average energy E vac = 0) if the quark 
masses are zero. In quantum chromodynamics (QCD) the vacuum is more complex with 
a finite mass of the quarks, which may be coupled to scalar Higgs particles, subject of 
extensive searches with the available and planned high energy experiments at Fermilab 
and CERN. 

To investigate hadron mass generation we start from a Lagrangian, in which two gluons 
couple to scalar fields, from which qq pairs are emitted. The possibility that two gluons 
coupled to + may produce colour singlet bound states has been mentioned already by 
Cornwall [2], but in context to the structure of QCD. Since the overlap of gluon fields is 
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non-local, such a model could serve also as an effective theory to investigate the strong 
scalar fields observed in hadron excitations [3] and scattering [I] [7], which may be difficult 
to extract from QCD. 

Assuming a scalar coupling of gluon fields of the form gg — > (qq) n we write the Lagrangian 
in the form 

Csi = - - A (F^ - G^) , (1) 

with D^ 1 being the covariant derivative = d^—ig s A^ and F^ u the Abelian field strength 
tensor F^ u = d^A u — d u A^, with being the gluon fields. couples two gluon fields 

with J w = + to gg-pairs with 

G" 1 = -g s A» [h{W lg ^) + b 2 (W lg V) 2 + & 3 (^i 9 *) 3 + •••] , (2) 

where (^Vi g ^/) stands for the creation of gg-pairs, which interact by gluon exchange. 
Corresponding Feynman diagrams are shown in the upper part of fig. 1. 

The structure of C s = ^G^G^ 1 implies a colour neutral coupling of two gluon fields. Hence, 
the symmetry is simply isospin SU(2) (two quarks with different charge and one gluon) 
without colour and (because of massless quarks as seen below) without flavour degree of 
freedom. By the coupling of two gluons to J 71 = + the Lagrangian has no chiral symmetry, 
leading naturally to the sequence of hadronic states as observed experimentally. 

A two-gluon field is produced only, if there is spacial overlap of two gluon fields. We write 
the radial part of C s by a matrix element 

$( ri - r 2 ) =a s < A a (n)| [&i(a}Vi a a ? ) + b 2 {a\V lg a q f + ...] 2 \A 2 (r 2 ) > , (3) 

where A^r) are radial gluon fields, and a q quark and antiquark creation operators, 
and V\ g an interaction dominated by 1-gluon exchange between these quarks^- For the 
following discussion only the first term in eq. (j3J) is needed, contributions due to (gg) 2 
contributions will be discussed briefly at the end of the paper. 

Due to the non-local 2-gluon field, given by a 2-gluon density p$(r), the recoiling local 
quark fields are also smeared out and a q Vi g a q can be written as a q-q potential V qq (r) 
described by folding a gg-density with the gluon-exchange interaction. The fact, that 
1 More-gluon exchange as well as spin-spin and spin-orbit effects have not been considered. 



2 



V qq {r) is a scalar potential but the created gg-pair has negative parity, requires a p-wave 
gg-density p p qq (f) = p qq (r) Y hm (6, (f>) leading to 

V qq (r) = J dr'ffjf') Y 1>m (6', 0') V lg {r - r') , (4) 

where V\ g (r) can act only within the density p$(r). Therefore the interaction has to be 
cut towards large radii and we use the following form 

Vi 9 (r) = {-a./r) e~ cr (5) 

with cut-off parameter c determined in a self-consistent fit of the 2-gluon density. Note, 
that eq. ([3]) indicates a departure from a purely relativistic description with a lifetime of 
the created system At. Causality is fulfilled, if At > (r — r')/c. The shape of the folding 
potential (jlj) has been calculated by multiplying the density and potential in momentum 
space and retransforming the product to r-space. 

The p-wave character of the gg-density gives rise to the constraint < r qq > = J dr rp qq {r) = 
and thus 

p q9 (r) = V3 (1 + ■ d/dr) p*{r) , (6) 

where (3 is determined from the condition < r qq >= 0. A consequence of eq. (j3J) is, that 
the radial dependence of V qq (r) should be the same as p$(r) 

Mr) = V qq (r) . (7) 

From the different relations between densities and potential in eqs. (HI) , ([6]) and ((7|) the 
two-gluon density is determined. Self-consistent solutions of eq. (E]) are obtained assuming 
a form 

p<f>(r) = p [exp{ — {r/a) K }} 2 with k ~ 1.5 (8) 

and interaction cut-off parameter c, which yields a mean square radius of the effective 
interaction between 20 and 80 % larger than < r| >. The self-consistency condition is 
rather strict (see the lower parts of fig. 1): using a pure exponential form (k—1) a very 
steep rise of p$(r) is obtained for r —>■ , but an almost negligible and flat potential, 
which cannot satisfy eq. (j7]). Also for a Gaussian form («=2) no consistency is obtained: 
normalised to the inner part of p<s>{r), the deduced potential falls off more rapidly towards 
larger radii than the density. Only by a density with k ~ 1.5 a satisfactory solution is 
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obtained. A self-consistent solution is also shown (in r- and momentum (Q)-space) in 
fig. 2 with a quantitative agreement for radii >0.1 fm and momenta <2.5 GeV/c . The 
localisation in space indicates, that a stationary (mesonic) system is formed. Interestingly, 
the emerging quarks can couple again to two gluons giving rise to stable glueballs as well. 

In the Fourier transformation to Q-space the process gg — ► qq in question is elastic 
and consequently the created gg-pair has no mass. However, if we take a finite mass of 
the created quarks of 1.4 GeV (such a mass has been assumed in potential models [8] 
for systems of similar size), the dashed line in the lower part of fig. 2 is obtained and 
no self-consistent solution is possible. Thus, our solutions require massless quarks and 
consequently the deduced hadronic systems can be related directly to the absolute vacuum 
of fluctuating gluon fields. 

For a quark-gluon system with finite mean square radius as shown in fig. [2] the corre- 
sponding binding potential can be obtained from a three-dimensional reduction of the 
Bethe-Salpeter equation in form of a (relativistic) Schrodinger equation 

V2/z$ 

where /i$ is the mass parameter, which is related to the mass m$ of the qq system by 
= \ m® — 8m with a correction 8m ^ only for very light systems. Further, ip^>(r) is 
the wave function, which is given for a system of uncorrelated gluons by ^^(r)! 2 = p$(r). 
The dependence of the potential as a function of ip&(r) is then given by 

= f- + -%^)-Vt + Vo • (10) 



( Ti l r d 2 2d 

-dr 2 r dr 



Vi(r)Wr) = Ei Mr) , (9) 



2/1$ ^ dr 2 r dr ' V ; *( r 
Inserting the form of the density in eq. (jHJ) yields explicitely 



VJr) 



h 2 



- 2 C-r~ 2 [K(-r-(K+i)]\+v . (ii) 

a a a J 



2/z$ 

Since the (2g — qq) system couples to the vacuum, V a is assumed to be zero. With 
the values of k and a from the self-consistent solutions ([8} and /i$ ~1.5 GeV (which is 
determined below also self-consistently) the deduced densities (shown in figs. [1] and [2]) yield 
a binding potential given in the upper part of fig. [3j This has the same form as the known 
confinement potential V con f = —a s /r+br deduced from potential models [8] and the lattice 
data of Bali et al. [9J . Further, the potential ffTTl) is the same for other solutions of smaller 
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or larger radii and can therefore be identified with the 'universal' confinement potential. It 
is important to note, that the self-induced confinement potential (fTTl) reproduces the 1/r + 
linear form without any assumption on its distance behavior; this is entirely a consequence 
of the deduced radial form (jBJ) of the two-gluon density. Interestingly, because of the rather 
simple confinement mechanism in our approach, a direct connection to strings (which have 
been related to the linear part of confinement potential) appears possible. 

To determine bound state energies of basic scalar and vector qq states (with J PC = 
++ and 1 ) binding in the self-induced confinement potential ( TTTT) but also in the q-q 
potential has to be considered, which (different from the confinement potential) depends 
strongly on the radius of the density. For smaller radii a deepening of this potential is 
observed, which leads to more strongly bound states. 

We define the mass of the system by the energy to balance binding (this is also the way 
hadron masses are observed). This yields 

Mi = —E qq + Ei , (12) 

where E qq and E{ are the binding energies in V qq (r) and V$(r), respectively (note that the 
eigenstates of V qq {r) have negative energy). 

We apply another constraint (generation of mass by binding), requiring that the mass M 
of the lowest bound state is equal to m$ (appearing in the mass parameter //$ = |m$ — 5m 
with 5m ~0). With this condition we find m$ ~ l/< r| > 15 for scalar states with a 
coupling constant a s decreasing for smaller systems. For the five systems in table 1 we 
find values of a s of about 1.3, 0.5, 0.4, 0.32 and 0.23. Together with other solutions we 
find a s (M Q ) « 1.5 af CD (Q) for M Q = Q up to large masses (with values of af CD (Q) 
from the systematics in ref. [10J). This yields evidence for asymptotic freedom also in our 
model. More details will be discussed elswhere. 

So far we have no constraint on the radial extent of the (2g — qq) system, giving rise 
to continuous mass spectra. Discretisation is provided by a vacuum potential sum rule 
^vac{f) = En $ n ( r )i requiring that the sum of two-gluon matrix elements ([3]) (where 
$ n (r) = ct™p$ n (r)) is equal to the total 1-gluon exchange force in the vacuum. For this 
we take the simplest form, demanding that the total cut-off by the finite size of the 
overlapping gluon fields shows the same 1/r dependence as the 1-gluon exchange force. 
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Table 1: Deduced masses (in GeV) of scalar and vector qq states in comparison with 
known ++ and 1 mesons [10]. Only the Is states in the q-q potential are given. The 
radii of the known states of mass M Q are fine-tuned to agree with experimental data. 



Solution 


(meson) 


M 


Mi 


M 2 


M e xp 


M exp 


Mf p 


1 ++ 


a 


0.60 


1.60±0.2 


2.30±0.3 


0.60±0.1 


2 ++ 


fo 


1.30 


2.15±0.2 


2.80±0.2 


1.30±0.2 






1 


Ijj 


0.78 


1.66±0.2 


2.28±0.3 


0.78 


1.65±0.02 




3 0++ 


fo 


2.30 


3.00±0.2 


3.55±0.2 


2.30±0.2 






1— 


$ 


1.00 


1.69 


2.27 


1.02 


1.68±0.02 




4 0++ 


not seen 


4.5±1.0 


M o +0.5 


M o +1.0 








1— 


J/ty 


3.10 


3.68 


4.19 


3.097 


3.686 


(4.160) 


5 0++ 


not seen 


22.0±8 


M o +0.4 


M o +0.7 








1— 


T 


9.46 


9.99 


10.36 


9.46 


10.023 


10.355 



From this we obtain 

rv 

^(r) = — = -E^W- (13) 

' n 

Five self-consistent solutions for scalar states below 50 GeV have been extracted with 
mean square radii < r\ > of about 0.49, 0.26, 0.12, 0.06, and 0.02 fm 2 , respectively, 
and masses M Q given in table 1. In the lower part of fig. 3 the different solutions $ n (r) 
for n=l-5 are shown together with their sum (solid line). This is compared to the sum 
rule (fl3|) with <5 S =0.35 (lower dot-dashed line). Good agreement is obtained, indicating 
that the sum rule is fulfilled. Only for radii <0.15 fm a deviation is observed, which is 
explained by contributions from solutions with masses significantly larger than 50 GeV. 

Between the scalar states discussed so far and corresponding vector states there is a direct 
relation: the p-wave qq density has to be replaced by the corresponding s-wave density. 
This leads to momentum distributions of the vector states shifted to smaller momenta 
(and corresponding smaller masses) than the scalar states. We obtain mean square radii 
< r 2 $ > of about 0.38, 0.24, 0.10, and 0.04 fm 2 and masses given in table 1, which are 
in good agreement with the masses of the strong 1 mesons (together with their radial 
excitations) of the (isoscalar) "flavour families" u, $, J/^f and T. The masses of the 
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known ++ states are also in general agreement with experiment: the lowest ++ state 
corresponds to the <r(600) meson, which has been clearly identified [H] as a broad meson 
resonance in J/^-decay. The next ++ states at 1.3 and 2.3 GeV may correspond to the 
scalar resonances / o (1300), see also ref. [12], and / O (2300), whereas the higher ++ states 
(not found so far) may not be observable in e + - e~ collision experiments. 



As compared to potential models with finite quark masses, as e.g. in ref. [8], we obtain 
significantly more states, bound states in the confinement and in the q-q potential. The 
solutions in table 1 correspond only to the Is levels in the q-q potential, in addition we have 
calculated Ns levels for N=2, 3, and 4. Most of these states, however, have a relatively 
small mass below 3 GeV. As the q-q potential is Coulomb like, it creates a continuum 
of Ns levels which ranges down in mass to the threshold region. This continuum should 
mix with the states in the confinement potential, giving rise to large scalar phase shifts at 
low energies, which are observed experimentally but so far not well understood in other 
models. 

Concerning masses above 3 GeV, solution 5 of table 1 yields additional ++ 2s and 3s 
states in the q-q potential at masses of about 12 and 8.8 GeV, respectively, whereas an 
extra 1 2s state is obtained (between the most likely \l/(3s) and \I/(4s) states at 4.160 
GeV and 4.415 GeV) at a mass of about 4.2 GeV. This state may be identified with the 
recently discovered X(4260), see ref. [10]. Corresponding excited states in the confinement 
potential ( fTTi) should be found at masses of 4.9, 5.3 and 5.5 GeV with uncertainties of 
0.2-0.3 GeV. 

By identifying the deduced qq states with known mesons we can check the overall con- 
sistency of our model by the observed widths of these states. For the lifetime we assume 
< At >=< r| > 1 / 2 jc and get for the heavier mesons in table 1 values which are small 
compared to the lifetime deduced from the relatively small experimental widths. Dif- 
ferently, for the lightest meson in table 1 we obtain a value of < At > of 2.4 10~ 24 s 
corresponding to a width of about 300 MeV, which is less than the width (~ 500-600 
MeV) extracted [TT] for a(600). However, as the Coulomb like q-q potential gives rise to 
a low energy continuum of Ns states, the width of the lowest bound state in table 1 has 
to be much smaller, indicating that also in this case our approach is valid. 

Finally, possible (qq) 2 contributions corresponding to the second term in eq. ([3]) are ad- 
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dressed. In the results shown in fig. 2 and the upper part of fig. 4 the folding potential 
deviates from p$(r) only at small radii. This difference could be filled by a small (qq) 2 
component. This is supported by Monte Carlo simulations [13] . in which good agreement 
between p$(r) and V qq {r) in eq. ([7]) is obtained in the entire r and Q region by the inclusion 
of a (qq) 2 contibution. 

In conclusion, a model has been presented, in which hadron masses are described as 
bound states of quarks with a very simple structure of the vacuum. This leads to a good 
description of the confinement potential and hadron masses. Other results, including a 
discussion of asymptotic freedom and the stability of baryons will be discussed later. 

We thank P. Decowski, M. Dillig (deceased), A. Kupsc and P. Raczka among many other 
collegues for fruitful discussions, valuable comments and the help in formal derivations. 
Special thanks to P. Zupranski for writing the Monte Carlo code, for numerous conversa- 
tions and clarificatons. 
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Figure 1: Upper part: Relevant Feynman diagrams 2g — > qq and 2g — > (qq) 2 (q denotes 
quark or antiquark). Lower parts: density p$(r) with < r 2 >=0.06 fm 2 (dot-dashed lines) 
and folding potential (solid lines) for k—1, 1.5 and 2, respectively. 
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Figure 2: Upper part: Self-consistent solution of eq. (J4]) for a scalar state with two-gluon 
density and q-q potential given by dot-dashed and solid lines, respectively. Lower part: 
Same as in the upper part but transformed to Q-space (multiplied by Q 2 ). The dashed 
line corresponds to a calculation assuming quark masses of 1.4 GeV. 
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Figure 3: Upper part: Confinement potential from lattice calculations [9] in comparison 
with the qq binding potential ffTTj) given by solid line. Lower part: q-q potentials for the 
different solutions in table 1 with their sum given by solid line compared to the vacuum 
sum rule ( 1T3|) . lower dot-dashed line. 
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